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The results of mathematical modeling of hot-cast molding of ceramic are presented. The computational data 
describe the changes of the aggregate state of thermoplastic slip as it cools in the molding cavity. The casting 
properties giving uniform properties of the thermoplastic slip during hot-cast molding of ceramic are obtained. 
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Hot casting of thermoplastic slip is currently the main 
method of molding complex parts as well as tubular and 
rod-shaped parts [1]. The essence of this fabrication method 
is to impart to liquid thermoplastic slip the configuration of 
the part by filling a molding cavity under pressure and then 
fixing the acquired shape by changing the aggregate state of 
the medium as it cools. 

The casting system consists of a multiphase medium 
where a thermoplastic binder with complex composition and 
coagulation properties is used as the dispersion phase [2, 3]. 
The rheological properties of the slip are determined by the 
ratio between the kinetically free and bound dispersion 
phase. Beryllium oxide is a promising material for a number 
of special applications in technology. However, its unique 
thermal conductivity at the molding stage gives casting sys¬ 
tems an elevated stiffness, which makes it difficult to control 
structure formation during slip motion [2]. 

The method of hot casting includes the following stages 
[1,2]: slip motion and heat transfer in the liquid state; slip 
motion and heat transfer taking account of solidification; cast 
motion and heat transfer in the solid-plastic state. At the 
stage when the molding cavity is filled and held under pres¬ 
sure it is most important to ensure the maximum destruction 
of the structure in order to obtain a uniform suspension [2]. 
During solidification it is necessary to attain minimum fric¬ 
tion against the wall of the molding cavity and maximum 
plasticity without repeated destruction of the newly formed 
cast structure. 

The hydrodynamics of thermoplastic slip during casting 
should be viewed as a physical process of deformation and 
flow. The slip flow retains its configuration after exiting the 
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feeder. Experiments have established that in the range of pos¬ 
sible casting rates the flow of thermoplastic slip in a casting 
mold is laminar [2]. Slip enters the casting mold at tempera¬ 
ture 75 - 85 °C and cools in the mold to 30 - 45°C, at which 
the cast can be extracted from the mold without being dis¬ 
torted [2]. 

The development and implementation of the method of 
hot casting have been largely empirical. A detailed analysis 
of the hydrodynamics and heat transfer will make it possible 
to approach on a more substantiated basis the technology of 
the slip casting regime. 

We shall examine the results of a mathematical model of 
the process of hot-cast molding of beryllium oxide ceramic. 

RHEOLOGICAL MODEL 
OF THERMOPLASTIC SLIP 

It is assumed that thermoplastic slip (high-viscosity sus¬ 
pension) is a two-phase disperse system, where the solid 
mineral phase is a beryllium oxide powder while the liquid 
phase is an organic binder [2]. The organic binder consists of 
three components: paraffin, bees wax, and oleic acid in the 
mass ratio (wt.%) 82: 15 : 3. Thermoplastic slip with 
thixotropic flow with shear rates in the range 0.005 - 
1200 sec -1 can be treated as a Shvedov-Bingham visco¬ 
plastic liquid [4]: 

ft/ 

+ 1- 1 > (!) 
or 

where t is the shear stress, Pa; x 0 is the maximum shear 
stress, Pa; jli is the plastic viscosity, Pa • sec; and, du/dr is the 
velocity gradient, sec “ '. 

The plastic viscosity jli and the maximum shear stress x 0 
of the slip depend on the temperature t. The experimental 
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Fig. 1. Variation of the thermoplastic slip density versus tempe¬ 
rature. 


data of [2] with relative mass content of the binder co = 0.10 
are described by the empirical relations 

p = 1.17 x 10 9 exp (- 0.280 + 49.4; 
x 0 = 2.61 x 10 10 exp(- 0.3390 + 22.34. 

The density of thermoplastic slip is determined by the 
beryllium oxide powder and binder concentrations 

p _ P BeO P cb _ 

(l-®)PcB + ®PBeO ’ 

where p Be0 is the density of beryllium oxide, g/cm 3 ; p b is the 
binder density; and, co is the relative mass content of the 
binder as a fraction of unity. 

The binder density depends on the temperature and is de¬ 
termined by the empirical relation 

P b (0 = 0.83 + 0.062cos [tt(0.035£- 5.15)]. 

The beryllium oxide density is p Be0 = 3.047 g/cm 3 . The 
binder density in the temperature range from 75 to 40°C va¬ 
ries within the limits p b = 0.773 - 0.9 g/cm 3 and, accord¬ 
ingly, on solidification the thermoplastic slip density p in¬ 
creases from 2.355 to 2.460 g/cm 3 (Fig. 1). 

Experimental data [2] show the empirical temperature 
dependences of the thermal conductivity X and specific heat 
c p of thermoplastic slip to be 

X = 1.6 + 4.8exp (- 0.017 1); 
c p = lOOOexp (t/ 290). 

In summary, the rheological properties of thermoplastic 
slip consisting of beryllium oxide are complicated functions 
of temperature and an aggregate change of the liquid suspen¬ 
sion into a solid-plastic body occurs during casting. 

MATHEMATICAL MODEL 

The subject of study is the motion and heat transfer of 
thermoplastic slip in an annular cavity (Fig. 2). A hot liquid 


217 


r 



Fig. 2. Flow in a coaxial cavity: r l9 r 2 ) inner and outer radii of the 
annular channel; r 3 ) radius of the outer wall of the channel; L x , L 2 , 
L 3 ) lengths of the hot, warm, and cold loops; q, t 2 , t 3 ) temperatures 
of the hot, warm, and cold loops, respectively. 

slip with initial temperature t 0 = 75°C flows into an annular 
cavity with thickness (r 2 - r x ) and length L. According to the 
experimental data [2] the maximum initial velocity does not 
exceed u 0 = 2 mm/min. The slip mass cools and solidifies as 
it moves, acquiring shape at the exit from the molding cavity. 
A distinguishing feature of beryllium oxide slip is high ther¬ 
mal conductivity, but the Prandtl number Pr = \xc p /X is much 
greater than 1 because of the viscosity of the thermoplastic 
slip. The wall of the molding cavity is cooled by water circu¬ 
lating in an annular case. The cooling zone is divided into 
three parts: the temperature t x = 73°C in the first part, 
t 2 = 59°C in the second, and t 3 = 45°C in the third. The total 
length of the cavity is L = 0.108 m and the lengths of the hot, 
warm, and cold parts are, respectively, L x = 0.022 m, 
L 2 = 0.045 m, and L 3 = 0.041 m. 

The intensity of heat transfer between the hot slip and 
cooling liquid depends on the flow regime, thermophysical 
characteristics of the slip, temperature of the cooling liquid, 
material of the tube wall, and geometric dimensions of the 
case. In the course of the motion the rheological properties 
change because the liquid slip solidifies. Heat is released on 
the transition surface of the aggregate state of the slip. The 
slip mass in stationary motion in an annular cavity is cooled 
by the water laving the tube on the outside. A characteristic 
feature of solidification in the cold zone of the molding ca¬ 
vity is that the temperature t t of the outer wall of the annular 
channel will be lower than the temperature t w of the inner 
wall with radius r ,; this imparts nonuniformity to the temper¬ 
ature profile as well as to the rheological properties of the ex¬ 
truded slip. Solidification starts on the outer wall, while the 
slip will be in a liquid state at the inner wall of the annular 
cavity. As a result, liquid slip may be added in order to com- 
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pensate the interior shrinkage at the stage of solidification of 
the cast in the cooled zone of the cavity. 

The problem is studied in a cylindrical coordinate system 
with axes z and r, directed parallel and radially to the flow, 
respectively. The slip motion is stationary in the direction of 
the casting velocity, and it can be studied using a system of 
equations of hydrodynamics, closed by the Bingham model 
of a non-Newtonian liquid. The thickness of the annular ca¬ 
vity is much smaller than the length, so that the narrow-chan¬ 
nel approximation can be used to study the motion and heat 
transfer of the slip [4]: 


du du 

p u -bpv — 

dz dr 


dp 

dz 



H- 

r dr 


r\x 

V 


du^ 
dr j 


+ 


1 d 

~ — (fT 0 ) + pg; 

r or 

dpu 1 dpv „ 
= 0 , 

dz r dr 


( 2 ) 

(3) 


dt 

puc — + pvc 

dz 


dt 

dr 


i a 


r dr 


Xr — 

TJ ^ 

+ EL + JLl 

r du' 

V dr) 

OS 

l dr J 



where z, r — axial and radial coordinates; u, v — compo¬ 
nents of the velocity vector; s — time; p, p, t, t 0 , c p , p, and 
X — pressure, density, temperature, maximum shear stress, 
specific heat, viscosity, and thermal conductivity of the slip. 

The heat of solidification H k of the beryllium oxide slip 
is found from the experimental data and equals H k = 
7800 J/kg [5]. The rheological properties of the slip are ex¬ 
pressed by empirical relations. 

The pressure gradient in the equation of motion is deter¬ 
mined from the conservation of the mass flow: 

2nj purdr = n(r 2 -r 2 )p 0 w 0 . (5) 

n 


The velocity and temperature distributions at the en¬ 
trance into the annular cavity are assumed to be constant 
over the cross section of the channel and, correspondingly, 
all thermophysical properties of the slip are constant: 

for z = 0: u = u 0 , v = 0, t = t Q . (6) 

Adiabaticity is assumed on the inner wall of the cavity: 

for z > 0: r = r { , dt/dr = 0. (7) 

The conditions of impermeability and slippage of the slip 
are imposed on the walls of the annular layer: 

for z > 0, r = r., v = 0,, / = 1,2. (8) 

Heat transfer on the outer wall occurs in accordance with 
the values of the temperature in the cooling loops of the ca¬ 
vity. Let t h t 2 , and t 3 be the temperature of the water in the 


hot, warm, and cold loops, respectively. The boundary condi¬ 
tions for the temperature on the outer wall are: 


dt 

for 0 < znl { , r = r 2 : -X — = k(t-t l ); 

dr 


(9) 


dt 

for / 1 < znl 2 , r = r 2 : —X — = k(t-t 2 ); (10) 

dr 


dt 

for l 2 < znl 3 ,r = r 2 : —X — = k(t—t 3 ). 

dr 


(ii) 


The coefficient of heat transfer k is calculated in the stan¬ 
dard manner: 


1 


1 1 
+-+ 
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kd 3 o i 2 d 2 CL 3 d 3 2X w \n(d 3 /d 2 ) 


where d x , d 2 — inner and outer diameters of the cavity, m; 
X w — thermal conductivity of the wall material, W/(m • K); 
a 2 , a 3 — inner and outer heat transfer coefficients, W/(m 2 • K). 

The coefficient of heat transfer a 2 is found from the ex¬ 
pression [6] 

q w l(dt / dr) w 


where q =-X 

± w 



heat-flux density at the wall, W/m 2 ; 


t m — average mass temperature; t w — cooling water tempe¬ 
rature. 

The average mass temperature of the slip in the annular 
cavity is determined in the standard manner [6]. 

Therefore, knowing the temperature and velocity distri¬ 
butions over the cross section of the cavity, the inner heat- 
transfer coefficient a 2 can be calculated. 

The outer heat-transfer coefficient a 3 depends on the 
properties of the cooling water and the dimensions of the 
case and is given by the Forchheimer relation: 



2^6 

d 3 In (4 d 4 / d 3 ) 


where d v d 4 — diameters of the annular case; X b — thermal 
conductivity of the cooling liquid. 

The system of Eqs. (2) - (11) is solved using dimension¬ 
less variables. The coordinates z, r are divided by the radius 
r,, the velocity components u , v by w 0 , the pressure p by the 
dynamic head, the temperature by x 0 , and the density, maxi¬ 
mum shear stress, specific heat, viscosity, and thermal con¬ 
ductivity by their values at the temperature x 0 . 

The system of equations (2) - (6) with the boundary con¬ 
ditions (7)-(11) is solved numerically. The computational 
region is divided into elementary cells with edges A z t , Ar •. 
The difference analogs of the equations of motion (2) and 
energy (4) are obtained using the second-order Crank-Ni- 
cholson scheme, and the difference analog of Eq. (3) is ob- 
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Fig. 3. Temperature (a) and velocity ( b ) distribu¬ 
tions in an annular cavity with t 0 = 75°C and 
Uq = 2 mm/min. 
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tained by a two-layer second-order scheme [7]. The pressure 
gradient is determined by the method of splitting [7] from the 
flow conservation condition (5). 

The problem of heat transfer in laminar flow of a Newto¬ 
nian liquid in a circular tube with a constant wall temperature 
was solved in order to verify the numerical method. The vari¬ 
ation of the Nusselt number along the tube was obtained in 
the calculations: it decreases monotonically, goes to a con¬ 
stant value Nu = 3.48, and is identical in accuracy to the ana¬ 
lytical solution [6]. 

DISCUSSION 

The calculations were performed for three values of the 
annular layer r 2 /r 1 = 1.32, 1.4, and 1.5. The entry slip velo¬ 
cities were u 0 = 1.0 mm/min and u 0 = 2.0 mm/min. The 
temperatures of the cooling water in the hot, warm, and cold 
loops were t t = 73, 59, and 45°C, respectively. 

Figure 3 shows the distributions of the velocity u and 
temperature t in different sections of the annular cavity 
r 2 /r x = 1.32 with casting velocity 2.0 mm/min. The inten¬ 
sity of the radial heat transfer, arising as a result of heat 
transfer at the outer wall, increases considerably the heat 
transfer from the hot slip to the cooling liquid. At the dis¬ 
tance z/r x = 0.75 it is evident that the gradient of the temper¬ 
ature head has a quality effect on the coefficient of heat emis¬ 
sion, just as the flow velocity gradient has on the slippage. 
Since the effect of the pressure gradient is the same over the 
cross section of the annular channel, the corresponding 
changes of the viscosity and maximum shear stress on the 
walls give rise to different values of the velocity on the cav¬ 
ity walls. The slip layers move relative to one another with 
different velocity, i.e., motion and heat transfer occur in the 
laminar regime. 

The temperature distribution is uniform over the cross 
section of the annular cavity (see Fig. 3 a). This is explained 
by the fact that because of the high thermal conductivity of 
beryllium oxide conduction plays the predominate role in 
heat transfer. Such a temperature distribution imparts unifor¬ 
mity to the rheological properties of the slip. The aggregate 
state of the slip changes at the horizontal section z!r x = 2.0. 


The heat released dissipates rapidly because of the high ther¬ 
mal conductivity of beryllium oxide. In the cold-loop zone 
the motion and heat transfer of the slip occur in the solid- 
plastic state. 

Thus, the uniform distributions of the temperature and 
rheological properties at all stages of molding preclude the 
formation of pits and porosity inside the ceramic as the slip 
solidifies. 

An increase of the thickness of the annular gap 
r 2 / r x = 1-5 results in a higher mass flow and heat content of 
the slip as well as in nonuniform temperature and velocity 
distributions over the cross section of the annular cavity 
(Fig. 4). The cooling rate of the slip decreases, and solidifica¬ 
tion starts later at the distance zlr x = 2.2 than in the case 
r 2 /r x = 1-32. The temperature nonuniformity over the cross 
section of the annular cavity results in nonuniformities in the 
distribution of the rheological properties of the slip. The 
phase transition starts at the outer wall. The heat of solidifi¬ 
cation of the slip increases even more the nonuniformity of 
the distributions of the temperature and, correspondingly, the 
rheological properties. This can cause the slip to shrink 
nonuniformly during solidification and internal shrinkage 
pits and porosity to form at the ceramic formation stage. 

Thus, nonuniformity of the temperature and, correspond¬ 
ingly, rheological properties of the slip can result in the for¬ 
mation of a defect during casting of ceramic parts. This is 
confirmed by the experimental data [2]. 

Figure 5 shows the variation of the density of thermo¬ 
plastic slip ( a ) and the heat flux density on the outer wall {b ) 
along an annular channel. The density variation in the hot 
loop is small (Fig. 5 a ). The density of the thermoplastic slip 
increases as the slip cools; the relative density changes are in 
the range 1-1.01 in the liquid state and 1-1.03 in the so¬ 
lid-plastic state. The computed data are in agreement with 
the experimental results [2, 3]. 

The variation of the heat flux density at the wall of the 
molding cavity is the determining factor in controlling the 
temperature regime of the method of hot casting of ceramic 
(Fig. 5b ). The heat flux density at the outer wall increases 
sharply in the initial part of the warm and cold loops and de¬ 
creases monotonically with decreasing slip temperature 
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Fig. 4. Temperature (a) and velocity (b ) distribu¬ 
tions in an annular cavity with t 0 = 75°C and 
Uq = 2 mm/min. 
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Fig. 5. Distribution of the slip ( a ) and heat flux ( b ) densities at the 
wall with t 0 = 75°C and u 0 = 2 mm/min. 

(Fig. 5b ). The zig-zag variation of the heat flux density is ex¬ 
plained by heat emission along the outer wall of the molding 
cavity due to heat removal in the hot, warm, and cold loops 
of the annular tube. 

The next regime was calculated for annular cavity size 
r 2 /r x = 1.4 and flow velocity u 0 = 1.0 mm/min. Compared 
with the two preceding regimes, in this case solidification 
starts at the distance z/r x = 1.9 and the phase transition front 
occupies a horizontal section. As the computational results 
show, this thermal regime also gives uniform rheological 
properties of the slip at all stages of molding, determined by 
the temperature field, and a ceramic part of good quality. 

CONCLUSIONS 

The following basic conclusions can be drawn from the 
results of this research. 

The mathematical model of ceramic molding by hot cast¬ 
ing includes a system of equations of motion of a non-New¬ 
tonian liquid, continuity, and energy taking account of heat 
release accompanying a change of the aggregate state as the 
thermoplastic slip cools in the molding cavity. 


The velocity, temperature, and rheological properties dis¬ 
tributions showing the internal structure of the slip in the 
molding cavity were obtained in the calculations. The com¬ 
putational data describe the entire stage of molding of a ce¬ 
ramic taking account of the change in the aggregate state and 
show that the uniformity of the temperature field over the 
cross section of the cavity with shrinkage of the thermoplas¬ 
tic slip plays the main role. 

The slip density increases as the aggregate state changes. 
The calculation explains the mechanism of the compensation 
of the volume changes during molding of the ceramic by hot 
casting and is in satisfactory agreement with experiment. 

The results of this study show that the mathematical 
model of the process is an effective tool for controlling the 
method of hot casting ceramic with the optimal conditions 
determined as a function of the design molds and technologi¬ 
cal parameters of the molding of thermoplastic slips. 
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